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We investigate the condensate wave function and elementary excitations of strongly interacting 
bosonic polar molecules in a harmonic trap, treating the scattering amplitude beyond the standard 
first Born approximation (FBA). By using an appropriate trial wave function in the variational 
method, effects of the leading order correction beyond the FBA have been investigated and shown 
to be significantly enhanced when the system is close to the phase boundary of collapse. How 
such leading order effect of going beyond the FBA can be observed in a realistic experiment is also 
discussed. 
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I. INTRODUCTION 



In recent years dipolar gases become a fast growing 
field of theoretical and experimental interests in the stud- 
ies of ultracold atoms and molecules. Among several 
dipolar systems, chromium (^^Cr) atoms have been suc- 
cessfully realized and studied to great extents [H-Q. The 
dipolar interaction effects for ®''Rb atoms as well as ^^K 
atoms are also observed in different groups 0, II] ■ Sev- 
eral polar molecule systems such as CO 0, ND3 p^ . 
RbCs [TTl, LiCs a nd CsCl have also been trapped, 
cooled, and studied More recently, a high phase- 

space density gas of polar ''°K*''Rb molecules have also 
been produced [l^ . These stimulate great interest in the 
studies of dipolar systems at low temperatures. In an 
earlier theoretical work within the first Born approxima- 
tion (FBA), Yi and You [11] proposed a pseudopotential 
to study the long-range dipolar interaction. Based on 
this approximation, which is justified only in the weak 
dipole moment limit, various theoretical studies of the 
excitations, collapses, instability, etc. of the dipolar BEC 
system have been carried out in these years \u\ . 

While Yi and You's pseudopotential within FBA is ap- 
propriate for a weakly interacting system, it can not be 
applied for a polar molecule system of strong interaction. 
Under a strong field, polar molecules can have a large 
electric dipole moment to make various shape resonances 
possible. Therefore one needs to include the interaction 
beyond the FBA in order to have a better understanding 
on the low energy behavior of these dipolar systems. Re- 
cently, one of us has developed an effective many-body 
theory for bosonic polar molecules in the strong inter- 
action regime that goes beyond the FBA [l^. It is in- 
teresting to investigate how the higher order interactions 
affect the properties of polar molecules. Following [Tsj , 
our current paper attempts to study the leading order 
effects beyond the FBA on the condensate profile and el- 
ementary excitations of bosonic polar molecules. It will 
be shown that the effect of the leading-order term is most 
prominent when the system is approaching the instability 



(collapsed state) . This is manifested in both the ground- 
state property and elementary excitations. In particular, 
at the intermediate trap aspect ratio, A = ujz/ojp = 5 ~ 6, 
the results obtained by the theories within FBA and go- 
ing beyond FBA have a drastic difference near the col- 
lapsed regime. This drastic change indeed allows one to 
extract quantitatively the effect beyond FBA in a real- 
istic experiment. In the large A (pancake) and small A 
(cigar) limits, on the other hand, the effects beyond the 
FBA becomes much reduced and may not be easily ob- 
served in the experiment. 

The paper is organized as follows. In Sec. II, we out- 
line the effective Hamiltonian for studying the low-energy 
many-body physics of polar molecules. In Sec. Ill, a trial 
wave function, called "modified Gaussian" is introduced 
and used to calculate the ground-state properties of polar 
molecules beyond the FBA. In Sec. IV, we investigate the 
elementary excitations (breathing modes) of the system, 
and provide an scheme to quantitatively extract the ef- 
fect beyond the FBA in an experiment. Sec. V is a brief 
conclusion. 



II. EFFECTIVE HAMILTONIAN 

The effective Hamiltonian, which describes low-energy 
many-body physics of bosonic polar molecules beyond 
the FBA, can be given in the following second quantiza- 
tion formalism (for details, see Ref. |18j): 
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where K!xt(pi — P2) is the Fourier transform of the ex- 
ternal trapping potential and VL is the system volume. 
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r(pi,p2) is the pseudo-potential responsible for the in- 
teraction vortex between polar molecules, and can be di- 
vided into the following three parts: 



r(p,p') = 



/a(p,p')-(2) 



M ^ ' M 

The first term in Eq. ([2]) is from the standard s-wave 
scattering and the second term is due to the usual FBA 
with Vrf(q) = 47r_D^(3cos^0g^ — l)/3 being the long-range 
dipolc-dipolc interaction. Here D is the electric dipole 
moment and 9q_, is the angle between the z-component 
wavevector and the total wavevector. The third term, 
/a, is the scattering amplitude including all the results 
deviated from the known FBA in the second term. In the 
general situation, the deviated term can be expanded in 
diiTerent angular momentum channels p^ . 



/a(Pi,P2) 
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where l^a\", ' = -V - tmm.'M is the differ- 

ence between a full scattering length and its FBA result. 
Yim{p) is the spherical harmonics defined by the orien- 
tation of momentum p. In Eq. ([3]), the sum, j, , has 
excluded the contribution from a pure s-wave channel 
{I = V = 0), which has been already included in the 
first term of Eq. In the limit of a weak external 

field, the electric dipole moment is also small and each 
term in Eq. ^ has the following orders of magnitude: 

as = 0{D''), Vd = OiD^), and Aa^f - 0{D^) 
Thus in the low field (i.e., small-D) limit, /a(p, p') term 
can be safely neglected within the FBA. However, when 
the external field is strong enough, the effect of /a(p, p') 
should be taken into account to go beyond the standard 
FBA. The actual values of Aa|™'' have to be calculated 
from a full scattering theory of polar molecules, as shown 
in Ref. [19-21]. 

In the present paper, we shall study how these higher 
order terms, Aa|™'', can affect the ground state and 
the elementary excitation property of a bosonic dipolar 
molecule gas beyond the first Born approximation level. 
According to the numerical calculation of these scattering 

amplitude in Refs. [l^[2l|, we find that Aal™"* becomes 
smaller for larger angular momentum, I or V . When near 
the first shape resonance, it is found that although all 
scattering channels will diverge, but the most dominant 
contributions are still from {1,1') — (0,0) (s-wave), and 
(0,2) = (2,0) channels. Divergences on the (2,2) chan- 
nels are almost invisible. When away from the shape 
resonance regime, the scattering amplitude of higher an- 
gular momentum channels {I > 0) are well described by 
the first Born approximation, if only the dipole moment 
is not too large. Therefore, from a practical point of view, 
we do not need to investigate all the possible values of 
Aai^\ but, on the other hand, we can concentrate on 
the effect of the lowest non-trivial terms, Aal^l ^ and 

set Aa'ij} = for aU {1,1') not equal to (0,2), (2,0), or 
(0,0). 



As mentioned above, this special case becomes relevant 
when the electric dipole moment is close to the first shape 
resonance peak, where the coupling between s-wave and 
d-wstve channels are significantly enhanced. The effect of 
shape resonance on the s-wave channel, Us, is also sig- 
nificant, but has been separated in Eq. When the 
dipole moment is even larger, it is reasonable to expect 
that more scattering channels will have scattering ampli- 
tude different from the first Born approximation result, 
i.e. Aal™"* ^ 0. However, since it is very difficult to in- 
vestigate a general properties of such strong interacting 
limit, here we will still concentrate on the effects beyond 
the FBA only through the channel, Aag'^2 Oj which will 
be a good approximation when the external field is near 
the regime of the first shape resonance [l^. Neverthe- 
less, the effect of next leading term, Aa^l, will be briefly 
studied in Sec. IV. 

To calculate the expectation value of the effective 
Hamiltonian Eq. E — (Hcs), one can replace Ok by 
a macroscopic condensate wave function, = (ok) = 
^ / (ir\E'(r) e^*'^ '" in ([T]) at zero temperature. We can 
then apply the variational method to obtain the ground- 
state energy of the system. If one uses the (normalized) 



simple Gaussian-type trial wave function {p 



*(r) 



exp(-pV2i?g - z^ml) 



1/2 



(4) 



with i?o and Rz the Gaussian radii of the condensate in 
the x-y plane and along the z axis respectively and as- 
sumes that the harmonic trapping potential is V^xt(r) = 
imw^p^ -I- imw^z^ with Wp and the trapping fre- 
quencies, variational energy of Eq. ([T]) becomes (see also 
Eq. (16) of [18]) 



E{Ro, Rz 

En 



E]i + Etrap + E\ 



Ant, 



where 



and 
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(7) 



^2(^)^8) 

^z 



correspond to the kinetic energy, trapping potential en- 
ergy, and interaction energy respectively. Note that the 
above results include only the leading higher-order term. 



a4% Here M^) 



Pdx) 



(l + ((fl.0/-R..)2-l)2:2)3/2 



and A = w^/wp is the trapping aspect ratio. All lengths, 
as, ad = MD'^/Ti', Aoq^j, Ro, and Rz, are scaled 
by the harmonic oscillator length, flosco = •\/7i/-/Vfa;p. 
i?o = Nh^ /ma'^g^ Q is the energy scale. Throughout this 
paper, we will use the same length and energy scales. 
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It is convenient to redefine: Cs = NasXd = Nad, 
and Co. 2 = —SV5NAa'j^2 such that the interaction term 
Eq. ([5]) can be rewritten as 



1 



C.|^ + ^(Cd + C0,2)A2(|^) 



.(9) 



In view of Eq. (jH) or ([8|), it poses a subtlety that both the 
dipolar interaction Cd and the higher-order interaction 
Co, 2 couple to the same function. This result is just an 
artifact of the Gaussian trial wave function and should 
not exist in a more general condensate wave function as 
stated in Ref. [Hi- 



III. BEYOND THE FBA: GROUND-STATE 
PROPERTY 
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In order to separate the contribution of the Aao''2 for 
the usual FBA result, in this paper we use the following 
"modified Gaussian" (MG) trial wave function 



*(r) = Cexp(- 



Rl 



2 



'ml) 



(10) 



where i?o, Rz, oq 
C = 2^-3/V[i^o^ 



and are variational parameters and 



Rz (Sag + 4aoa3 + Saf ^ Soq + Aaz + 4)] 
is the normalization factor. Function (jlOp can be viewed 
as the simple Gaussian (SG) multiplied (modified) by a 
parabolic function. In Ref. [22| , a similar MG trial func- 
tion (SG function multiplied by a hypercosine function) 
was first introduced and used to study one- and two- 
component ultracold BEG systems. When the MG trial 
wave function, Eq. (ITUl) . is used to calculate the energy 
fimctional, E = E{Rq, i?^, oq, a^), ground-state energy of 
the system can be obtained by minimizing E. That is, it 
requires that dE/dPi = are satisfied for all four vari- 
ational parameter /3i = Rq, Rz,ao,az. Comparison be- 
tween the results due to a SG and a MG trial wave func- 
tion for one- and two-component ultracold BEG systems 
without a dipolar interaction can be found in Ref. 



22| 



To justify the appropriateness of MG trial wave func- 
tion pO)) in the variational studies of strong dipolar BEG 
systems, we refer to the following. Recent numerical 
studies on dipolar BEG systems have revealed struc- 
tured ground-state density profiles in them [2^, [23| • Of 
most interest, density profile can exhibit a double-peak 
structure, so-called "biconcave condensates", for dipo- 
lar BECs. In a recent work of Jiang and Su •29], it has 
been shown that the structured double-peak ground state 
density profile can be successfully reproduced in the vari- 
ational approach based on a similar trial wave function, 
analogous to (|10l) . Thus it is believed that, at least in the 
qualitative manner, MG trial wave function should be a 
good one for studying the higher-order effect beyond the 
FBA in strongly interacting dipolar BEG systems. It is 
hoped that a more accurate numerical approach, which 
is difficult but is underway, can be completed soon. 



FIG. 1. (Color online) In frame (a)-(c), the condensate aspect 
ratio, 7 (see text for the definition), is plotted for MG against 
the value of dimensionless Cd (Co, 2 set to zero) or Co, 2 (Cd set 
to zero) for three trapping aspect ratio A = Uz/uup = 0.35, 5 
and 10. Frame (d) shows the phase boundary separating the 
unstable (collapsed) state from the stable state. In all frames, 
Cs = is taken. 



When the MG (|10p is used as the trial wave function, 
the three terms of the variational energy become 



Rl 
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ARlRl 
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Rq + Rl 



RlRl 
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- 8ao 
2ao - 



- 

'daz 



h 12a3 + 3aj 
2a^ + 15a? 



h 4aoa2) 
- SaoCz)], 



and 



-Eint = Cs^s (flO, Oz: ^0, Rz) + QAd (oq, ttzjRo, Rz) 

+ Co, 2^0,2 (ao,a^,i?o,^2) , (H) 

where B = {Auz + Suq + 3al + 8ao + 4aoaz 4-4). Ag, Ad, 
and Ao^2 are lengthy functions of ao,az,Ro, and Rz whose 
explicit forms are given in Appendix [X] One can simply 
check that when ag = Uz = 0, B ^ A and -Ek and i?trap 
reduce to those for the SG case [see Eqs. ^ and ([7])]. E'int 
will also reduce to that for the SG case (see Appendix El) . 

As seen clearly in Eint in Eq. ([TTl), Crf^d and Co,2^o,2 
correspond to the contributions of dipolar interaction and 
higher order AaQ*^2 interaction to the ground-state energy. 
Thus the regime at which the two terms become most 
distinct is also the regime to see the effect beyond the 
FBA most clearly. In this section, we will investigate such 
beyond-FBA effect from the condensate aspect ratio, 7, 
which is defined as 



7 



(12) 
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where 



|*(r 



|2 ^2 



dr 



|*(r)|^a;2^r = / |*(r)|^ y^dr. 



(13) 



In Fig. [Tla)-(c), based on the MG trial wave func- 
tion, we show results of the ground-state aspect ratio, 7, 
against the dipolar interaction strength, C.j_, for 2 — 0, 
and against <^o,2 by setting Q — for three different val- 
ues of the trapping aspect ratio, A. Note that, although 
such kind of interaction is not realistic in the experiment, 
but it provides a direct evidence to distinguish the effects 
of dipolar interaction in the FBA and beyond the FBA. 
In view of Fig.[lja) and (c), for the A = 0.35 and 10 cases, 
the effects of and Co, 2 are almost the same within the 
parameter regime we calculate, because the condensate 
aspect ratio is mostly determined by the trapping aspect 
ratio directly, i.e. the interaction effect is negligible. In 
contrast, for the case of A = 5 [see Fig. (Hb)], the effects 
of Cd and Co, 2 are quite distinct: the system becomes col- 
lapsed more easily for the case of a pure FBA interaction, 
Q (solid line), while it becomes less easily collapsed for 
the other case (dashed line) . In Fig. [IJd) , we show the 
critical value of C for these two cases as a function of the 
trapping aspect ratio, A. One sees that system becomes 
always stable when A is larger than 5.3 for the bare FBA 
result, while it requires 6.5 for a pure Aag''2, the leading 
order effect beyond the FBA. Our result indicates that 
Co, 2 acts more attractively than Cd in the intermediate A 
regime. 

It is useful to note that regarding where the system 
collapses, there might be some quantitative difference 
between the variational and numerical results. Under 
the condition of same A, real system could collapse more 
easily than what the variational theory predicts in the 
smaller dipolar interaction regime [23l |26| . In addition, 
for large A (the system is of pancake shape), the system 
could also be unstable at large (d 22b i2M- Nevertheless, 
we emphasize that the qualitative features predicted by 
the variational method should still be reliable for large Cd 
so long as A is not too large. The latter is exactly what 
is studied in the current context. 



IV. BEYOND THE FBA: ELEMENTARY 
EXCITATIONS 

This section devotes to the elementary excitations of 
a dipolar condensate beyond the first Born approxima- 
tion. As sketched in Fig.[2l two types of collective modes, 
namely the in-phase and out-of-phase breathing modes, 
are considered here. They are of particular importance 
because they are the lowest two excitation modes in a 
trapped condensate [H, . In the paper of Goral and 
Santos [13], these two modes were also called breathing 
and quadrupole modes. For simplicity of presentation, 
we will show the results only for the lower-energy one. 



t 

in-phase mode 



1^ 



out-of-phase mode 



FIG. 2. Scliematic plot of the in-phase and out-of-phase 
breathing modes of an anisotropic dipolar gas. In Ref. [T?!, 
these two kinds of modes are called breathing and quadrupole 
modes. 



although both of them are simultaneously obtained in 
the same method. 

In this paper, variational method is applied to study 
the collective excitation of a system. It corresponds 
to solving the stationary point of the action S = 
J dtL, where the Lagrangian L — T — E with T ~ 
J dv{ih/2)[^*{r)d-^{v)/dt--^{r)d-^*{v)/dt] and E is 
obtained by the sum of the terms in Eq. (fTTj) . For 
the study of breathing modes, the MG trial function in 
Eq. (jlOp will be generalized to include dynamical vari- 
ables as follows: 

P^(l + .0(t)+i4(t)) ^2(l + e^(t) + ,e'^(t)) 

* (r, i) = C (t) e L 

, ao(l+eo(t))/52 , a,{l + e,{t))z'^ 



7?2 



(14) 



Here, in a cylindrically symmetric trap, £i and e[ {i — 
0, z) correspond to the fluctuations of local amplitude 
and local phase of the dipole cloud associated with the p 
and z directions. As mentioned before, the values of i?Oj 
flo, and az are determined by minimizing the energy 
functional. After some lengthy derivations, we obtain the 
dispersions for the breathing modes: 

= {/1/4 + /2/5 ± [(/1/4 - f2hf + ^flhhY'^}/^, 

(15) 

where ± correspond to the in-phase or out-of-phase 
modes and 



, 1 a^s 

/l — 77 



, 1 a^s 

/2 — 77- 



fi = B(4+ 16ao + 4a2 



24a^ 



/3- 
f 3a' 



d^E 

dR^dRa 

1 . ^-l 



(16) 



/s ^ 2B(4 + 8ao + 12a^ + 8a^ + 15a^ 



To determine which one corresponds to the in-phase or 
out-of-phase mode in (I15|) requires solving explicitly the 
time dependence of the dynamical variables. 
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FIG. 3. (Color online) Frame (a) and (b): Lower breathing 
mode frequencies are plotted as a function of A for three com- 
binations, (Cd,Co,2) = (10, 10), (20,0), and (40,-20). Here 
C,s = Q for (a) and = 2 for (b). Frame (c) and (d): Parallel 
to frame (a) and (b), 7 ratios are plotted as a function of A. 
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FIG. 4. (Color online) (a) Lower breathing mode and (b) 7 
ratio of the dipolar system are plotted as a function of (2.2- 
In both frames, A = 6 and two combinations of (Cd,Co,2) = 
(10, 10) and (40, —20) are chosen. Here we set (s = 0. 



In the following, a case study is given first, which will 
lead to an explicit idea how to extract quantitatively the 
^0,2 term in real experiments. In ^^Cr atom dipolar BEC, 
it was measured that D'^M/h^ ~ 24 A. If atom number 
is about N = and the harmonic oscillator length 
is about 1/im, C,d will be about 20. For polar molecule 
systems, could be 100 times larger than that of ^^Cr 
atom dipolar BEC [27]. Nevertheless, here we study the 
case of Q + Co, 2 = 20 to which three combinations of 
iCd, Co,2) = (10,10), (20, 0), and (40, -20) are considered. 

In Fig. ^a) and (b), the lower breathing modes for 
these three combinations are presented and compared as 
a function of trapping aspect ratio A. In the current con- 
text, the lower breathing mode oscillates out-of-phase. 
This is mainly due to the fact that s-wave interaction 
is relatively small and the dipolar interaction dominates. 
Since dipolar interaction is anisotropic: attractive along 
the z direction while repulsive in the xy plane, conse- 
quently out-of-phase mode has lower energy than that of 
the in-phase mode. Fig. M^a) corresponds to the case of 
Cs = 0. In view of Fig. |3Ia), two important features are 
revealed: (i) all three curves merge in the large A limit 
(i.e. a pancake like trapping potential) and (ii) the curves 
deviate each other significantly when the system is close 
to the phase boundary of collapse. For instances when 
A = 5.9, w = 1.338^0 for the case (Cd,Co,2) = (20,0), 
while w = 1.044a;o for the case {(d, (0,2) = (40, -20). The 
relative frequency difference, Aw = \uji — lu2\/uj2, can be 
21% large. Aw will increase even more significantly when 
the system is approaching the collapsed regime. Similar 
behaviors are also found in Fig. [2Kb), where a finite value 
of s-wave scattering length is included {(s = 2). 

In Fig. [Hfc) and (d), we show the condensate aspect 
ratio, 7, as a function of A. In view of Fig. EJc) with 
Cs = 0, one finds that all three 7 curves also merge in 
the large A limit. Besides, when A = 5.9 close to the 
collapsed regime, 7 = 0.353 and 0.392 respectively for the 
case of (Cd,Co,2) = (20,0) and (40,-20). The relative 7 



difference, A7 = [71 — 72I/72, is about 11%. In Fig.|3Kd) 
with (s = 2, when A = 3.8, 7 ratio is 0.528 and 0.593 for 
the case of (Cd, (0,2) = (20, 0) and (40, -20). This gives a 
relatively smaller A7 — 2%. Similar to the Aw case, A7 
will increase significantly when the system is even more 
close to the collapsed regime. 

It should be emphasized that in the case that (g is large 
or the sum of + Co. 2 is small, the system will tend to 
stabilize over a large span of A. This means that whatever 
combinations of Cd and Co, 2 under Cd + Co, 2 = const will 
roughly lead to the same curve. Thus, Aw or A7 studied 
above, will always be small, i.e. the effect of the effects 
beyond the FBA cannot be observed easily. 

In the following, we propose how to extract quantita- 
tively the value of Co, 2 by measuring the lower breathing 
mode frequency or the condensate aspect ratio 7 of the 
system. Firstly, we assume that the value of Cs is rela- 
tively small such that the behaviors of breathing mode 
frequency and 7 ratio resemble those in Fig. [2Ia)-(d). 
Suppose that one does not know the value of any one of 
Cs, Cdj and Co, 2- One can first perform the measurement 
of the lower breathing mode frequency and/or 7 ratio at 
large A case, say A = 10. As shown in Fig. [Sla)-(d), all 
three curves merge at large A limit. This means that one 
can unambiguously determine the values of both Cs and 
the sum of Cd and Co, 2 simply by carrying out a theoret- 
ical fitting to the experimental data. The next task is to 
separate the values of Cd and Co, 2- For this purpose, one 
can redo the experiment and decrease A to approach the 
collapsed regime. Since at this regime, the fitting will be 
quite sensitive to both values of Cd and Co, 2, one can then 
unambiguously determine the value of Co, 2 by comparing 
with our theoretical calculation result. 

Now we provide a short discussion on the effect of 
the next order effect beyond the FBA by including finite 

value of Aa^l- Same as before, we define C2,2 = A^Aa2''2 
for the convenience. In Fig. 3] we show the lower breath- 
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ing mode frequency and the aspect ratio, 7, as a function 
of C2,2 for a fixed trapping aspect ratio, A = 6. Two com- 
binations of (Cd,Co,2) = (10,10) and (40,-20) are cho- 
sen and it is assumed that = 0. As ^2,2 increases, it 
is found that the system becomes more cigar-hke shape 
(i.e. 7 becomes larger), while it becomes more pancake- 
like shape when (^2,2 is negative. Positive C2,2 will lead 
the system to be more attractive along the z direction 
and hence more easily collapsed. However its effect is 
relatively weaker compared to that of the Co, 2 term dis- 
cussed before. 



V. CONCLUSION 

This paper attempts to study the ground state and 
elementary excitations of a strongly interacting dipolar 
bosonic gas based on a theory going beyond the first 
Born approximation (FBA). By using an appropriate 
trial wave function in the variational method, the lead- 
ing higher-order corrections to the FBA are studied in 
details, in particular for the condensate aspect ratio and 
the elementary excitation mode frequency. Of most in- 
terest, it is found that the effect of the corrections is most 
important when the system is close to phase boundary 
of collapse. Through several case studies with parame- 
ters pertaining to ^^Cr atom dipolar BEG, it is believed 
that the higher-order effect going beyond FBA should be 
even more significant and highly observable in the strong 
dipolar molecule system. An idea for extracting quanti- 
tatively such leading-order effect beyond the FBA in real 
experiments is provided. Finally, to shed more light on 
understanding the higher-order effect, a more accurate 
numerical calculation is in demand. 
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Appendix A: Forms of Functions A 

In this Appendix, we give explicit forms of functions 
Ad, and Ao,2 appearing in Eq. (|TT]) . Function Ag is 
analytic and given by 
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^ "32""^ ' 2 

Functions Ad and ^0,2 are represented by an integral in 
the cylindrical coordinate, (kr,kip,kz). More precisely 
Ad is given by 



An 



dkzj 2'Kkrdkr rjd 




(A2) 



with 



Vd 



{Riki+Rtki)/2 (3fc2/(^2 ^ ^2) _ 1) 



127r2 (4 Oo 02 + 8 Oo + 8 + 4 + 4 02 + 3 a2)2 
X [4 Oo 02 + 8 Oo + 8 Oo + 4 + 4 02 + 3 02 
- (2 Oo + 2 + 3 02) a2i?2fc2 - (2 + 4 Oo + 02) aoRlk"^ 



alRtkt alR^k^ oo a, R^Rlk^X 



4 4 
while Ao,2 is given by 



with 



00 CX3 

^0,2 = J dkzj 2iTkrdkr 770,2 

-00 



^fc^i?.gfc;)/2(3fc2y(fc2^fc2)_^) 



(A3) 



(A4) 



-{li 



Vo 2 — 

127r2(4 Oo 02 + 8 Oo + 8 + 4 + 4 02 + 3 a2) 

X [CrzaottzRlRlk^kl + DrttlR^k^ + DzalR^k'^ 

+ Co + CraoRlk"^ + CzttzRlkl] . (A5) 

Here Co, Cr, Cz, Crz, Dr, and Dz are given as follows: 

Co = 40 o^ + (44 02 + 96) og + (33 + 96 O2 + 112) og 



+ (|o? + 60al 



Cf — — 



— o^ + 18 o^ -f 34 al + 32 02 + 16, 
/21 

( 3 2 2 2 

I — O2 + 11 Oo 02 + 13 02 + 14 ciqUz + 32 Oo 02 
02 + 40 o^ + 48 Oo + 16 + 16 ojj) , 

I 1 c; 



15 

— Oo o^ + 18 0q02 + 28 O2 



3 ^ 2 15 

+ 36oo az + 20 ag + 40 o^ + 40 oo + 16) , 
Cr,z = 00^+402 + 200 02 +8 + 800 + 40^^ , 

Dr = + 202 + 00 02 + 4 + 4oo + 2o^^ , 

+ 4 + 4oo + 2o^j. (A6) 



/3 2 

Dz ^ [ —az + 2 02 + Oo a 



When Oo = ^ 
become equal. 



= O2 = 0, the two integrals for Ad and A^ 



-0,2 
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